#######################################################
# 
# --- Command log file ---
# PLAXIS 2D 2025.1
# Tutorial - 14
# Stability of dam under rapid drawdown
# Note: Use the PLAXIS Commands runner to execute this
# 
#######################################################
_initializerectangular SoilContour -130 -30 130 30
_setproperties "Title" "Stability of dam under rapid drawdown" "Company" "" "Comments" "" "UnitForce" "kN" "UnitLength" "m" "UnitTime" "day" "UnitTemperature" "K" "UnitEnergy" "kJ" "UnitPower" "kW" "WaterWeight" 10 "ReferenceTemperature" 293.15 "LiquidSpecificHeatCapacity" 4181.3 "LiquidThermalConductivity" 0.0006 "LiquidLatentHeat" 334000 "LiquidThermalExpansion" 0.00021 "LiquidTemperature" 293.15 "IceSpecificHeatCapacity" 2108 "IceThermalConductivity" 0.00222 "IceThermalExpansion" 5E-5 "VapourSpecificHeatCapacity" 1930 "VapourThermalConductivity" 2.5E-5 "VapourSpecificGasConstant" 461.5 "UseTemperatureDepWaterPropsTable" False "ModelType" "Plane strain" "ElementType" "15-Noded"
_borehole 0
_soillayer 0
_setsoillayerlevel Borehole_1 1 -30
_soilmat "Identification" "Core" "SoilModel" "Mohr-Coulomb" "DrainageType" "Undrained B"  "Colour" 15262369 "gammaUnsat" 16 "gammaSat" 18 "ERef" 1500 "nu" 0.35 "EInc" 300 "sURef" 5 "sUInc" 3.0 "verticalRef" 30 "GroundwaterClassificationType" "Hypres"  "GroundwaterSoilClassStandard" "Very fine" "GwUseDefaults" False "PermHorizontalPrimary" 0.0001 "PermVertical" 0.0001
_soilmat "Identification" "Fill" "SoilModel" "Mohr-Coulomb" "Colour" 10676870 "gammaUnsat" 16 "gammaSat" 20 "ERef" 20000 "nu" 0.33 "cRef" 5 "phi" 31 "psi" 1 "GroundwaterClassificationType" "Hypres" "GroundwaterSoilClassStandard" "Coarse" "GwUseDefaults" False "PermHorizontalPrimary" 1 "PermVertical" 1 
_soilmat "Identification" "Sub-soil" "SoilModel" "Mohr-Coulomb" "Colour" 10283244 "gammaUnsat" 17 "gammaSat" 21 "ERef" 50000 "nu" 0.3 "cRef" 1 "phi" 35 "psi" 5 "GroundwaterClassificationType"  "Hypres" "GroundwaterSoilClassStandard" "Coarse" "GwUseDefaults" False "PermHorizontalPrimary" 0.01 "PermVertical" 0.01 
_setmaterial Soillayer_1 Subsoil
_gotostructures
_polygon (-80 0) (92.5 0.0) (2.5 30.0) (-2.5 30.0)
_cutpoly (-10 0) (-2.5 30)
_cutpoly (10 0) (2.5 30)
_setmaterial (Soil_2) Core
_setmaterial (Soil_3 Soil_4) Fill
_mergeequivalents Geometry
_gotomesh
_mesh 0.04
_gotowater
_activate Polygons InitialPhase
_set InitialPhase.Identification "High reservoir"
_set InitialPhase.DeformCalcType "Gravity loading"
_set InitialPhase.PorePresCalcType "Steady state groundwater flow"
_waterlevel (-132 25) (-10 25) (93 -10) (132 -10)
_setglobalwaterlevel UserWaterLevel_1 InitialPhase
_rename UserWaterLevel_1 "FullReservoir_Steady"
_phase InitialPhase
_set Model.CurrentPhase Phase_1 
_set Phase_1.Identification "Rapid drawdown"
_set Phase_1.DeformCalcType "Fully coupled flow-deformation"
_set Phase_1.TimeInterval 5
_duplicate FullReservoir_Steady
_rename UserWaterLevel_1 "FullReservoir_Rapid"
_headfunction
_set HeadFunction_1.Name "Rapid"
_set Rapid.Signal "Linear"
_set Rapid.Time 5
_set Rapid.DeltaHead -20
_setglobalwaterlevel FullReservoir_Rapid Phase_1
_set WaterSegment_4.TimeDependency "Time dependent"
_set WaterSegment_4.HeadFunction Rapid
_phase InitialPhase
_set Model.CurrentPhase Phase_2 
_set Phase_2.Identification "Slow drawdown"
_set Phase_2.DeformCalcType "Fully coupled flow-deformation"
_set Phase_2.TimeInterval 50
_set Deform_2.ResetDisplacementsToZero True
_duplicate FullReservoir_Steady
_rename UserWaterLevel_1 "FullReservoir_Slow"
_headfunction
_set HeadFunction_1.Name "Slow"
_set Slow.Signal "Linear"
_set Slow.Time 50
_set Slow.DeltaHead -20
_setglobalwaterlevel FullReservoir_Slow Phase_2
_set WaterSegment_7.TimeDependency "Time dependent"
_set WaterSegment_7.HeadFunction Slow
_set Model.CurrentPhase InitialPhase 
_phase InitialPhase
_set Model.CurrentPhase Phase_3 
_set Phase_3.Identification "Low level"
_set Deform_3.IgnoreUndrainedBehaviour True
_waterlevel (-132 5) (-60 5) (93 -10) (132 -10)
_rename UserWaterLevel_1 "LowLevel_Steady"
_setglobalwaterlevel LowLevel_Steady Phase_3
_set Model.CurrentPhase InitialPhase 
_phase InitialPhase
_set Model.CurrentPhase Phase_4 
_set Phase_4.Identification "High reservoir - Safety"
_set Phase_4.DeformCalcType "Safety"
_set Deform_4.UseDefaultIterationParams False
_set Deform_4.MaxSteps 30
_set Model.CurrentPhase Phase_1 
_phase Phase_1
_set Model.CurrentPhase Phase_5 
_set Phase_5.Identification "Rapid drawdown - Safety"
_set Phase_5.DeformCalcType "Safety"
_set Deform_5.ResetDisplacementsToZero True
_set Deform_5.UseDefaultIterationParams False
_set Deform_5.MaxSteps 50
_set Model.CurrentPhase Phase_2 
_phase Phase_2
_set Model.CurrentPhase Phase_6 
_set Phase_6.Identification "Slow drawdown - Safety"
_set Phase_6.DeformCalcType "Safety"
_set Deform_6.ResetDisplacementsToZero True
_set Deform_6.UseDefaultIterationParams False
_set Deform_6.MaxSteps 50
_set Model.CurrentPhase Phase_3 
_phase Phase_3
_set Model.CurrentPhase Phase_7 
_set Phase_7.Identification "Low level - Safety"
_set Phase_7.DeformCalcType "Safety"
_set Deform_7.ResetDisplacementsToZero True
_set Deform_7.UseDefaultIterationParams False
_set Deform_7.MaxSteps 50
_gotostages
_activate Polygons InitialPhase
_selectmeshpoints
/output addcurvepoint "Node" Soil_2_1 -2.5 30
/output addcurvepoint "Node" Soil_3_1 -80 0
/output update
_set (InitialPhase.IgnoreSuction Phase_3.IgnoreSuction) False